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Abstract 

We study the impact of leading soft-collinear effects on threshold and joint-resummed 
calculations for the production of prompt photons via parton fragmentation, complement- 
ing a previous study for direct production. We assess these effects for both fixed-target 
and collider kinematics. We find them to be small, but noticeable and comparable to the 
direct case. 



Deceased. 



1 Introduction 



The perturbative QCD description of prompt production at hadron colliders can involve 
sizeable corrections from soft and collinear parton emission. In particular, the presence of 
a threshold at fixed pt induces large logarithmic corrections [H[2]. Expressed in terms of a 
(Melllin) moment variable N with the property that as N — > oo the kinematics approaches 
the threshold limit, these corrections take the form (L = IniV), 
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where the aij depend in general on the process. Such large logarithmic corrections can be 
organized and controlled through all-order resummation, for threshold [HHH] and joint [9)413] 
resummation. 

In recent years, in the context of threshold resummation, other large classes of terms have 
been brought under all-order control, such as large constants (' £ 7r 2 terms") |14H16j and a series 
of so-called "soft-collinear" terms 
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In a previous study [TTj we performed the resummation of the leading terms (j = 2i — 1) for 
direct production of prompt photons. In this paper we extend this study to the leading soft- 
collinear terms (which we shall also refer to as "In N/N terms") in the complimentary case 
of prompt photon production by fragmentation. In this case more subprocesses contribute, 
and their color states must be accounted for. These issues occur as well in threshold resum- 
mation without soft-collinear effects for both mechanisms of prompt photon production [18] . 
Subleading soft-collinear terms have been studied recently in other contexts in Refs. [19)431]. 

The paper is organized as follows. In section 2 we review briefly the threshold and joint 
resummed prompt photon pr distribution, and discuss the inclusion of soft-collinear effects. 
In section 3 we assess the numerical impact of these corrections, and we conclude in section 
4. In three appendices we clarify various technical points. 



2 Resummed transverse momentum distributions 

We consider the inclusive transverse momentum distribution of prompt photons produced at 
fixed pt in hadron-hadron collisions at center of mass (cm) energy 

Ha(pa) + h B (p B ) -»• l{p c ) + X , (3) 

where Ha b refers to the two incoming hadrons and X to the unobserved part of the final 
state. The lowest order QCD processes producing the prompt photon directly at partonic cm 
energy yfs are 

q(Pa) + q(Pb) -> 7(Pc) + gipd) ^ 

g(pa) + q(pb) -> jipc) + q(pd) , 

where in the second reaction q stands for both quark and anti-quark. The mimimum invariant 
mass s required for the final state is kp\. It is convenient to express the distance above 
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threshold by the variable 1 — x T , where x\ = Ap T /S. At the parton level this becomes 
1 - x\ = 1 - 4p T /s. 

Apart from the partonic sub processes that directly produce the photon, there are contri- 
butions from 2 — > 2 parton scattering 

a(p a ) + b(p b ) -)■ c{p c ) + d(p d ) , (5) 

where the photon is produced by fragmentation of final state parton c. In this paper, given the 
accuracy to which we work, this will be either a quark or anti-quark. The fragmentation com- 
ponent also contributes at 0(aa s ), as does the direct component, though it is subdominant 
in the sense that the fragmentation function behaves as 1/N [5]. This is in part because, for 
pp and pN collisions, the fragmentation component proceeds via valence quark scattering, as 
opposed to direct component which involves either a gluon or a sea quark. Morever, threshold 
resummation can substantially enhance this fragmentation component [18]. Here we consider 
the contribution of the fragmentation component to the threshold and joint-resummed cross 
section for prompt photons at fixed px when also the leading soft-collinear effects are included. 
The resummed cross section consists of two parts 

3 j (resum) 3 , (direct) 3 , (frag) 

p T aa AB +x p T ao AB ^ +x "t AB-**t+X 

~A = A + A ( 6 

apT dpT dpT 

where the two terms correspond to the subprocesses dH and ©, respectively. 

The expression for the joint- and threshold-resummed pt distribution of the direct compo- 
nent of the prompt photon hadroproduction cross section was derived in |10j (expressed in a 
somewhat different form in p2]). To be able to compare the expression for the fragmentation 
component below in Eq. (|15p with that for the direct component, we quote the latter result 
again here, but explain its structure only briefly. To next-to-leading logarithmic accurary, 
with leading In N/N effects included, it reads 



p T da 



AB^-y+X _ Pt 



dpT - ^E/J/g?/^e.-„,-,Q,l) 

x I dx T (x T ) N ^-^k^L C ab ^ d (a s (fi),x T ) (- S ) 

Jo Jl-xl V4|pt-Qt/2|V 



'^ a/A (Q, b, N) C b/B (Q, b, N) exp [^ T (JV, b, fx, Q) + E? T (N, b, /x, Q)] 



x exp 



F d (N,Q,»)+gW d (\) . (7) 



Joint resummation resums threshold and recoil logarithms together in terms of the variable N, 
the Mellin conjugate to x T , and the impact parameter b, the Fourier conjugate to Qt [9lfl0]. 
The latter is the recoil transverse momentum of the underlying hard scattering process, over 
which in the top line the integral is taken. The variable fl is a cut-off on this transverse 
momentum. The hard scale Q is in the present case equal to 2pT, and [i is the renormalization 
scale. 

The second line contains a Mellin transform over the partonic scaling variable x T in the 
recoil frame (indicated by the tilde), the Born amplitude M ab , the N- and 6-independent hard 
virtual corrections C a b_). 7 rf(a s (//), x T ) and a kinematic factor linking recoil (through Qt) and 
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threshold (through N) effects. Finally, the last two lines contain the Sudakov exponentials 
from initial and final partons as well as soft wide angle radiation in combined (N, b) space. 
They also feature the coefficients C a /^(Q, b, N) and C b / B (Q, b, N) which contain the evolution 
matrix for evolution of the parton distribution functions fh/ii{N, (jlf) from scale fj,p to scale 
Q/x and the parton distribution functions. 

The initial state perturbative exponent is given by 

Ea T (N,b,Q,fi) = *h®(J3) + h<P(J3,Q,ii) , (8) 
a a {n) 

with 

P = b a s (n) In (x) ■ (9) 

The functions fti ' 1 ^ are listed in the Appendix [Aj The function x(6Q, AT) defines the A and 
b dependent minimum scale of soft gluons to be included. As in [T7] we choose it to be 

^■ N) = i + T^jN- (10) 

where 

N = Ne~ tE , b = bQe lE /2, (11) 
7s being the Euler constant. The functions C a /^(Q,b, N) and C b / B (Q,b, N) are given by 

C h/H (Q,b,N) = J2 £ hg(N,Q/ X ,»F) f g/H {N^ F ) . (12) 
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where the matrix £ is the evolution matrix which implements evolution from scale \ip to scale 
Q/x- in so doing [32] one includes the leading In N/N effects due to initial state radiation in 
the evolution kernel [T7j. Note that, as a consequence, the expressions for ha and given 
here are slightly different from the versions in Ref. [TO]. To the accuracy we work, the fip 
dependence essentially cancels in Eq. (fT2|) . 
The final state exponent reads [33] 

F d (N, Q, (j,) = -^fi° ] (A) + fi 1] (A, Q, M ) + f' k (\, a 8 ) (13) 

where 

A = 6 a s (/i 2 ) In A (14) 

The functions fj^\x) , /^(A) and f' k (X), as well as fi^(A) in Eq. (jTj) , the exponent due to 
wide angle soft radiation, are all listed in Appendix lAl 

If one wishes to sum only threshold-enhanced logarithms, an easy modification of Eq. (JT|) 
|10yi7| suffices. One merely neglects the recoil term Qx in the kinematic factor in the third 
line of Eq. ([7j), upon which the Qx integral sets b to zero. 

Having reviewed the resummed cross section for the direct component of prompt photon 
production let us now turn to the fragmentation component. The joint-resummed expression 
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for the fragmentation component was derived in Refs. 



. It is given by 
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4|pt-Qt/2| 2 

C a6 ^ C(J (a s (yLi),Xr 



Af+1 



x exp 



ab—>cd 



N,b, 



■ /. 'I 



(15) 



Here D^/ C (2(N — 1) + 3,fi F ) is the fragmentation function which expresses the probability 
of parton c fragmenting into a photon. To include the soft-collinear effects we shall now 
make modifications similar to those for the direct component. There is however an impor- 
tant difference associated with the fragmentation function, which we shall discuss further 
below. Soft-collinear effects due to radiation from initial state partons a and b, can again 
be included by replacing the parton distribution functions by the functions Ci/^Q, b, N) and 
Cj/B(Q,b, N), and then of course modifying the associated exponent functions accordingly. 
To include the soft-collinear effects due to radiation from parton d, we shall modify the final 
state exponent as in Eq. (]13p . Finally, for parton c, associated with the fragmentation func- 
tion, one could envisage two approaches. In the first, one can include the soft-collinear effects 
directly via a modified exponent function, in analogy to the treatment for parton d. This is 
the approach we shall study numerically in section [3l In a second approach one could attempt 
to include them instead through evolution, in analogy to partons a and b. We discuss this 
approach, which is problematic, further in appendix ICl 

The modification required for the first approach is derived as follows [X7 l l35 | [36] . Starting 
with the usual integral form for the initial state exponent in threshold resummation 



E% (N, Q) = 2 / dz 
'0 



1 z n-i_ 1 rQ(i-*)dfi 



1 - z 



— A(a s (ii) 



one replaces, for a final state fragmenting (anti)quark, 



JV-1 



JV-1 



1 



JV-1 



1 - z 1-z 
As a consequence, carrying out the integral, one finds 

E? T (N, Q, fj,, fx F ) = £ c PT ' cik (iV, Q, fj,, ftp) + E^ SC (N, Q) , 



(16) 

(17) 
(18) 



where the first term contains the leading and next-to-leading threshold logarithms, and the 
second term the leading soft-collinear effects. The result is 

1 



E 



PT.eik 



(N,Q,/i,n F ) 



«s(m) 



,(0) 



(A) + g«(A) 



(19) 



where A = boa s lnN. The functions q^°'^(X) are defined in appendix [Al For the soft-collinear 
effects we find 

E^> sc (N,Q) = -^e W ' A 



2vr6 



b a s 



ln(l - 2A) 



(20) 
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We now arrive at the expression 



„3 ^( fra s) 
i>T aa AB^-y+X 

dp T 



abc JL 

x /l^r J ******** &-Qt) 



s 



N+l 



(2vr) 2 J ^ \4p' T 2 

x C l/A (Q, b, N) C j/B (Q, b, iV)4™ } ( N ~ 1, & ) 



(21) 



Here, S^^™^' (iV, 6) is the resummed cross section for the partonic process ab — > cd in com- 
bined N, b space and reads 



= exp [E* T (N, b, Q, fjt) + E? T (N, b, Q, + E? T (N, Q, /i, /i F ) + F d (iV, Q, »)] 
x Tr J H(Q,fi)Pexp 



Vt/N j I 

PT 



xPexp 



PT ' N dpi' 



PT 



ft 



(22) 



Here E^ T (N, b, Q, //), E? T (N,b,Q,n) correspond to Eq. flgj), F d {N,Q,fi) to Eq. ([13} and 
E^ T (N,Q, /j,, up) to Eq. (fT8j) . The exponents inside the trace are associated with wide angle 
soft radiation and its form is discussed below. 

The expression for the resummed exponent differs crucially from the corresponding ex- 
pression for the direct component in that there are a larger number of color structures that 
can connect the external partons in ([5]) as opposed to direct production. This requires a 
set of corresponding coefficient functions and soft anomalous dimension matrices, which may 
mix under soft emissions. These have already been computed for the threshold-resummation 
studies of p]|J[34|- In appendix [B] we recall the derivation of these factors. As a well-known 
result of this mixing of different color structures the radiative factor for wide angle soft ra- 
diation takes the form of a matrix in the space of allowed tensors that connect the color 
representations of partons a, b, c, d into a singlet. The trace in Eq. (|22p is taken in that color 
tensor space. S is a soft gluon function that represents non-collinear soft gluon emission, 
while H is the hard scattering function describing the short distance hard scattering. Both 
H and S are matrices in color tensor space. At lowest order, Sli = Tr[c^c/] where and cj 
are color tensors [39]. The soft anomalous dimension matrix T$ represents the evolution of 
the soft function from scale pt/N to px- The symbols P and P denote path ordering in the 
space of color tensors. 

Clearly, in a color basis in which the soft anomalous dimension T$ is diagonal, the path- 
ordered exponentials of matrices in Eq. (122p reduce to a sum of simple exponentials. Expres- 
sions for these soft anomalous dimension matrices have been given in Ref. [39] in terms of the 
mandelstam invariants s,t,u associated with the 2 — > 2 kinematics of reactions ([4]) and ([5]). 
In the threshold limit x T — > 1 one can approximate these invariants in the soft anomalous 
dimension matrices by 



t 



u 



-2p 2 T . 



(23) 
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We will illustrate this diagonalization procedure for the case of qq — > qq in Appendix[Bj After 
the diagonalization procedure is carried out, the resummed exponent for a given partonic 
channel is given by 



(iV-l,6)=C a ^ cd exp 



E^ T (N, b, Q,fx) + Ef T (N, b, Q, M ) 



+ E? T {N, Q, 11, n F ) + F d {N, Q, n) 



E G ^exp(r^W 
L / 



(24) 



The sum runs over all possible color configurations / with G ab ^ cd representing a weight for 

each color configuration such that Y2i ^ib-tcd = The anomalous dimensions r^^*^(iV) are 
given by 



r 



I,(int) 



ab-^cd 
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1 Z N-1 
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Di,ab->cd (a s ((l- z) Q 



(25) 



The NLL expansion of ^^^(N) is given by 



r 



I,(int) 
abated 



(AO 



D 



(i) 

Lab-^cd 



2-Kbn 



ln(l-2A) + 0(a s (a s In AT) 



(26) 



The coefficients D ab _^ cd , the color weights G 1 , one loop hard coefficients C ab _^ cd and Born 

cross sections in N- space sj^^ have been given in the Appendix of Ref. [34], for each of the 
partonic subprocesses in ([5]). Eqs. (|2T1) and (p4^) are the expressions we used for the results 
in section [3l 

Having presented the expressions for joint and threshold resummed production of prompt 
photons including soft-collinear effects, for both direct and fragmentation component, we next 
examine the latter numerically 



3 Numerical studies 

In this section we study numerically the inclusion of the In N/N terms for the case of prompt 
photon production for two kinematic conditions: those of pp collisions at the Tevatron at 
\/~S = 1.96 TeV [40p41| . and those of the pN collisions in the E706 |42| fixed target experiment 
with -Ebeam = 530 GeV, corresponding to = 31.5 GeV. As in our previous study [TT] . 
we stress that our aim is primarily to assess the effect of such terms in relevant kinematic 
conditions, rather than perform a comprehensive comparison with data. Our assessments 
therefore mainly consist of comparing the same calculation with and without In N/N terms, 
and compare this with the difference between LL and NLL accuracy. 

Our default choices for various input parameters are as follows. We use joint resum- 
mation unless specified otherwise. We use the GRV parton density set [43], correspond- 
ing to a s {Mz) = 0.114, with the evolution code of Ref. |44] . changing flavor number at 
p, = m c (1.4 GeV) and nib (4.5 GeV). We choose the factorization and renormalization scale 
equal to pr, and the non-perturbative parameter #np in Eq. ([33]) equal to lGeV 2 . For the 
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Figure 1: Relative contributions vs. px of direct and fragmentation photons for LL (a), NLL 
(b), NLL + In N/N (c). The left pane shows the results for E706 kinematics and the right 
pane for Tevatron kinematics. 



parameter x we use the expression in Eq. (152H . following [45], with r\ = 1/4. For our joint- 
resummed results, we chose for Tevatron (E706) kinematics the cut-off \x in Eq. (fit)]) equal to 
15 (5) GeV. Regarding logarithmic accuracy, and unless stated otherwise, we refer to LL when 
using only h ( a\ /f } and q(°\ and C (ab ^ = C (ab ^ = 1 (see Appendix [A]); we refer to 

NLL when also including f%, cf^ and the virtual corrections discussed below Eq. (I4ip . 
When including the soft-collinear In N/N terms we use the full NLO anomalous dimension, 
evolved from scale ^ to Q/x in Eq. (|46l) . 

In Fig. [T] we show the relative contributions of the direct and fragmentation components 
for both kinematic conditions to the total result 

direct fragmentation , , 

2 (27) 

direct + fragmentation ' direct + fragmentation ' 

for LL, NLL, as well as NLL with soft-collinear effects. For E706 kinematics we observe that, 
as pt increases, the contribution from the direct component becomes more dominant, but 
not such that the fragmentation component becomes wholly negligible. We note in particular 
that in the case of NLL with In N/N the relative direct contribution decreases slightly over 
the range of pt values. 

In the case of Tevatron kinematics, for values of pt < 150 GeV the relative direct contibution 
is larger than the relative fragmentation contribution. The contribution of soft-collinear effects 
is mostly larger than the difference between LL and NLL, the latter difference being quite 
small for joint resummation. In Fig. [2] we show the fragmentation component by itself for 
both kinematics, and for our three approximations. In the case of E706 kinematics we note 
that the NLL curve is decreased slightly relative to the LL curve. Furthermore, inclusion 
of the soft-collinear contribution leads to a further small decrease. In the case of Tevatron 
kinematics we see that the difference between the curves is even smaller, due to the smaller 
value of a s at larger px values. Inclusion of the In N/N terms lowers the result with respect 
to the NLL curve for px above 20 GeV by a few percent, but their contribution is negligible 
for pt values above 110 GeV. In Fig. [3] we show a comparison between joint- and threshold 
resummation for both E706 and Tevatron kinematics. The curves are all relative to the joint 
NLL result without soft-collinear terms. While for E706 and joint resummation the inclusion 
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Figure 2: The LL (a), NLL (b) and NLL+lniV/JV (c) calculations as well as the ratios 
NLL/LL, NLL+ln iV/iV/NLL vs. p T . Left pane: E706, right pane: Tevatron. 
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Figure 3: Comparison of joint resummation and threshold resummation effects, ratios to 
joint resummation NLL without In N/N. The three curves show joint NLL with In N/N (a), 
threshold NLL without In N/N (b) and threshold NLL with \nN/N (c). Left pane: E706, 
right pane: Tevatron. 
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Figure 4: Same as the ratio plots in Fig. [2] but for threshold resummation. Left pane: E706, 
right pane: Tevatron. 



of the In N/N terms leads to only a small decrease of the NLL curve , we see that the effects 
are much more appreciable for threshold resummation. In particular we note that for px 
above 4 GeV the inclusion of \nN/N terms leads to a notable decrease of the NLL curve. 
This is very similar to what we found for the direct component in Ref. |17| . For Tevatron 
kinematics the inclusion of In N/N terms for joint resummation is only appreciable for quite 
low values of pr- For threshold resummation the situation is again notably different, with 
the inclusion of the soft-collinear terms leading to an appreciable increase relative to the pure 
NLL result, again in similarity to the direct component |17j . Finally, Fig. [J] is equivalent 
to the ratios shown in Fig. [2] but now for threshold resummation only. The soft-collinear 
terms change from increasing to decreasing the pt distribution for E706 kinematics, with the 
opposite pattern occurring for Tevatron kinematics. 

The behavior in Figs. [3] and 2] is quite similar to the direct case in our previous paper, which 
is perhaps not surprising since the resummation effects (LL, NLL and In N/N) are introduced 
by means of process independent exponentials and evolution of parton distribution functions. 
The extra contributions from the fragmentation parton and the additional color structures 
do not seem to change these patterns noticeably. 

4 Conclusions 

We have examined the effects of including terms of the form 

^1 ^ 1 3 AT 

^EV^T- (28) 

3 

in the production of a hard photon through fragmentation for px distributions at both col- 
lider and fixed target kinematics, at leading accuracy (j = i), both for joint and threshold 
resummation. This is the complement to an earlier study carried for the direct production 
component |17j . 

To the extent that leading terms of the form (|28|) arise from initial state radiation effects, 
we used the method of Refs. [321145] to include them. Those arising from final state emission 
we included both by extending the resummation exponents for the fragmenting parton and 



10 



the unobserved parton to leading In N/N accuracy. As for the direct component, we found the 
combined In N/N terms to be comparable to NLL corrections, and dependent on kinematics, 
either enhancing or suppressing. The contribution from In N/N terms is noticable in the case 
of joint resummation although it remains small. This seems to imply that the corrections 
introduced by recoil effects tend to overshadow those of the soft-collinear terms. In the case 
of pure threshold resummation the soft-collinear effects play a more appreciable role. Overall 
we found a behavior comparable to the direct case studied previously. 
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A Exponents 



Here we list the exponents used in section [2j The initial state exponents for the LL and NLL 
case without inclusion of the In N/N— terms are given by 



^i 0) (A,/3) 
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(30) 



while for the case when evolving the parton distribution functions down to Q/x (Eq. (|H|)) 
they are 



/4 0) (/3) 
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where j3 = bQa s (n) lnx, and 



a(i) - a - ±r 



2 V 

with C q = C F and C 9 = Ca- Also we have 



/ 67 vr 2 \ 10 



(33) 



sW = -|c J r, j B« = -7t6 . (34) 



In these equations 



, IIC A -AT R N F nC 2 A -10C A T R N F -6C F T R N F 

b ° = IfcF ' 6l = 24^ • (35) 

where Tr = 1/2. 

The functions q^iX) in Eq. {EJ are obtained by setting /3 = A in Eq. ([29]) and (l30l) . 
The final state exponents (|13p involve the functions 



^ 0) = " " 2A) ln(1 " 2A) " 2(1 " A) ln(1 " A)] (36) 
7 » 1) = " 4^ 1 [ ln ( 1 - 2A )- 21n ( 1 - A ) + ^ ln2 ( 1 - 2A )- ln2 ( 1 - A )] 

+ 2€ ln(1 - A) -2^ [21n(1 - A) - ln(1 - 2A)] 

+ ^-[21n(l-A)-ln(l-2A)]ln^ (37) 

with A = boa s (fj>)ln N . Note that f3 only differs from A for joint resummation; for threshold 
resummation f3 = A. The final term in Eq. (|13p reads 

A (1) / A \ 

^=2k^{-^) lHl - 2X) - Hl - x)l ' (38) 

/ A \ 

4 = -f- exp [ln(l - 2A) - ln(l - A)] . (39) 

y 2-KbQ \ a s b J 

The wide-angle soft radiation exponents in Eq. ([7]) are 

9aiW = -%ln(l - 2A)ln2, <?« (A) = -%ln(l - 2A)ln2 (40) 



These expressions are obtained by expanding the perturbative functions A a (a s ), Bd(a s ) and 
Dab^d-f in powers of a s 

^K) = -4 1} + (-) 2 4 2) +0(^). (41) 

The explicit forms of 

are shown in [U[5]. The C^ ab ^ cd ^ are calculated by expanding 
the resummed cross section (Eq. (12ip ) to O(af) and matching to the fixed order NLO result 
in [36]. We note that there is no factorization scale dependence in and the coefficient 



functions in Eq. ([7]) because of complete evolution from scale fi F to Q/x m Eqs. ([T4"]) . ([2"3~ 
and CG]). 
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B Diagonization of soft anomalous dimensions 

In this appendix we discuss the diagonalization procedure for the anomalous dimension ma- 
trices of Ref . [6] to arrive at the expressions of G I c j J _ ¥cd and D\j_^, of Ref. [3l] , which we also 
use here. In what follows, we illustrate this procedure for the partonic sub process qq — > qq 
in the limit (f23|) . 

When one changes to a new color basis in which T$ is diagonal the resummed cross section 
takes a simpler form. The diagonalization procedure has been discussed in detail in Ref. [6j. 
Upon diagonalization the trace in Eq. (I22p reduces to a simple sum of exponentials as in 
Eq. (|24p . The soft anomalous dimension matrix for the process qq —> qq expressed in the 
£-channel singlet-octet basis is given by 



T s > = — 

TT 



2C F U -±(T + U) 2U 
fy-U 2C F T 



(42) 



where 



U = In ( — ) + wr, 



In 



ITT . 



(43) 



This results in a slight difference from the initial state exponent in Ref. [18] which we are 
using. The difference for the partonic sub process ab — > cd is then given by 



^ (C« + C h - C c + C d ) In (1-2A) 
which we add to the diagonal terms of anomalous dimension matrix. This leads to 

2 



(44) 



l s , = — 

7T 



9 



The eigenvalues of this matrix (|45[) are Ai = and A2 
vectors are, respectively 

-3 



1 



In 2. (45) 
= —2, and the corresponding eigen 

(46) 



Now, we change the basis from i-channel singlet-octet basis to a new basis in which Ts is 
diagonal by using the matrix of eigenvectors 



TV 



1 1 



(47) 



Thus we obtain 

Tr \h{ Pt , fi)?exp 



PT / N dn' 



rV f {a s (//)) 



PT 



xS(ajm]p 



Tr \ W l H(p T , [i)TlTV~ x exp 



— T s > la 8 l/i 

p T M v v 



Pt/N dfJ , 

~r r s 



PT 



fj,' 



TVRT 1 S exp 



dfi 



I L S 



n 



(48) 
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For the present case, 




2 



(49) 



H and S are the Mellin moments of 2 x 2 hard and soft matrices in color tensor space [6J. 
Substituting for them, one finds that the trace finally reduces to a sum of exponents 



o"i exp 



-Ailn2 
2nb 



ln(l - 2A) 



+ <r 2 exp 



-A 2 In 2 
2vr6 



ln(l - 2A) 



(50) 



where &\ and &2 are H^S'^ and H' 2i S' i2 respectively, H' and S' being the hard and soft matrices 
in the new basis. Identifying ai and -22— with Gl„,„„ and G 2 .„„ and and ~^ 2 , ln2 

with D qq _^. qq and D qq ^ qq , one obtains the relevant term in the exponent in Eq. ([24"]) . 



C Soft-collinear effects in photon fragmentation function 



In this appendix we discuss the possibility of including the leading soft-collinear effects in a 
manner exactly analogous to the initial state in |17j . i.e. Q/N. The present case is however 
special because the evolution equation for the photon fragmentation function is inhomoge- 
neous. 

To leading order, the non-singlet evolution equation for the photon fragmentation function 
Dy/ C (N, fip), where c is the parton that fragments into a photon, reads, in moment space 

where to 0(1/N), for 3 active flavors 



k<®(N) 



4 

3~/V ! 



p(°)(JV) =C F ( -2]nN 



J_ 3 

iV + 2 



(52) 



The inhomogeneous term arises because the photon, in contrast to hadrons, has a pointlike in- 
teraction with a quark. Correspondingly, the solution to Eq. (I5ip is the sum of a homogeneous 
(hadronic) and inhomogeneous (pointlike) part 



D\N, l ?) = Dl &d {N^ 2 )+Dl(N^ 2 ) 
To leading logarithmic accuracy the homogeneous solution is 



L 



p(°)(iV) 



where 



L 



— = 1 - 6 a s (/i )ln — 



(53) 



(54) 



(55) 



To see what evolution from Q to Q/N corresponds to, we substitute fio = Q and fi = Q/N 
and find 



L 



p(°)(JV) 



1 + 



exp 



2& a s ln(iV))^ ln(i ^ 



_ 1) 
2N it 



a s CF 



7T 



21n 2 (iV ) + ^-^ln ( iV) 



(56) 
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which indeed sums the leading soft-collinear effects. 

We now discuss the point-like (inhomogeneous) solution to Eq. (|5ip . To LL accuracy it 
reads 



1 - L 



1 p(0) (jy) 



1 



o 



^ (0) W 8vr 2 6 

2-Kbo 



k {0) {N) 



(57) 



Notice that D\{N, /j.^) = 0, and that it is indeed proportional to 1 /N via k^(N), see Eq. ([52 
For Ho = Q and ^ = Q/iV the factor in square brackets can be rewritten as follows 



1 rl-TT--P (0) W 

1 — L 2,ri, o v y 



p(°)(jy) 

&0 



;i + 2b a s ln(iV)) 



p(°)(JV) 



(58) 



We can rewrite the contents of square brackets as 



a s C F 



7T 
7T 



21n 2 (iV) + 



ln(iV) 



N 



■ ln(iV) 



(l + 2& Q a s m(iV)) 



2C F ln 2 (N) + C F 



ln(N) 
N 



-C F ]n(N) + 2-Kbo ln(iV) 



(59) 



The expression in Eq. (157D then becomes 



P" I .Y 9l 
pi \ ' N 2 



2C F ln 2 (N) + C F 



ln(AQ 
AT 



C F ln(iV) + 2vr6 ln(iV) 



a 



p(°)(JV) 



8tt 2 & 



fc(°>(JV) L 2 -"o (60) 



Carrying out the division inside the brackets for large N yields 



ln(N)k {0) (N)L 



,(0) 



(N) 



Taken together the photonfragmentation function takes the form 



pi { N 91 

^had I iV ' jy2 



1 — L l 2^b 

3vr N 



(61) 



(62) 



Let us now perform a rough assessment of the validity of this approach to include soft- 
collinear effects. As the hadronic component of the photon fragmentation function has the 
property that one can sum the leading In N/N effects through evolution, we factor it out 
of the full solution in (|62p . and estimate the large iV behavior of the remainder based on a 
reasonable assumption of the non-perturbative TV-dependence. We parametrize the hadronic 
part (stemming from the light quark flavors) in this limit as in Ref. |47j by -Dhad ~ /cix~°' 3 (l — 
x) 2 + A; 2 Vi(1.703 - x). It is now straightforward to show, after a Mellin transform, that the 
term for the point-like part contributes effectively at large N by a multiplicative factor 



l + a s (C\nN+9-lnN) 



(63) 



to the hadronic component for some constants C, C . This affects the leading \nNjN term in 
|, and would therefore seem to make this approach problematic. 
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